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ABSTRACT 

Oxygen abundances in the spiral galaxies expected to be richest in oxygen are esti- 
mated. The new abundance determinations are based on the recently discovered ff 
- relation between auroral and nebular oxygen line fluxes in high-metallicity Hii re- 
gions. We find that the maximum gas-phase oxygen abundance in the central regions 
of spiral galaxies is 12-|-log(0/H) ~ 8.75. This value is significantly lower (by a factor 
> 5) than the previously accepted value. The central oxygen abundance in the Milky 
Way is similar to that in other large spirals. 
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1 INTRODUCTION 

The oxygen abundance in spiral galaxies has been the sub- 
ject of much discussion for a long time. Which spiral galaxy 
is the oxygen-richest one? And how high is the oxygen abun- 
dance in the oxygen-richest galaxies? Due to the presence of 
radial abundance gradients in the disks of spiral galaxies, the 
maximum oxygen abundance occurs at their centers. The no- 
tation (0/H)c = 12+log(0/H)fl=o will be used thereafter 
to denote the central oxygen abundance. 

The chemical composition of various samples of spi- 
ral galaxie s has been discussed in a number of pa- 
pers ' " 
[ 1994 

l2002f) . According to those articles, the oxygen-richest 
lalaxies are: NGC 5194 (M 51) with (0/H)c = 9.54 



(e.E. IVila-Cost 


as & Edmundsl 119921 IZaritskv et al 


iGarnett et al.l 


Il997l 


Ivan Zee et alJ Il99gl: iGarnett 



JVila- Costas fc Edmund: 



Jll992^ ■ NGC 3351 with (0/H)c 

9.41 JZaritskv et alj 119*93^ NGC 3184 with (0/H)c 
9.5 (Ivan Zee et alJ '199'^ NGC 6744 with (0/H)c = 
9.64 llGarnB^^^^T^ 97l. Different versions of the one- 
dimen sional empirical method, proposed first bv lPaeel et alJ 
^23) a quarter of a century ago, have been use d for oxy- 
gen abundan ce deter minations in those papers. IPilvugii] 
J200d l2001allbl l2003l) has shown that the oxygen abun- 
dances in galaxies determined with the one-dimensional 
empirical calibrations are significantly overestimated at 
the high-metallicity end ( 12 -|- log O/H > 8.25). This 
for two reasons. First, the then-existing calibrating points 
at the high-metallicity end were very few and not reli- 
able. Second, the physical conditions in If 11 regions cannot 
be taken into account a c curat ely in one-dimensional cal- 
ibrations. IPilvugin et alj i2004l) have used instead a two- 
dimensional parametric empirical calibration to derive oxy- 



gen abundances for a sample of spiral galaxies which includes 
NGC 3184, NGC 3351, NGC 5194, and NGC 6744. They 
found generally lower oxygen abundances with (0/H)c ~ 
9.0. 

Here we consider anew the problem of the maximum 
oxygen abundance in spiral galaxies by attempting to derive 
more accurate oxygen abundances. These can be derived via 
the classical Tc - method, Tc being the electron temperature 
of the H 11 region. Measurements of the auroral lines, such as 
[OIII]A4363, are necessary to determine Te. Unfortunately, 
they are very faint and often drop below the detectabi lity 
level i n the spectra of high-metallicity H 11 regions. Pilvu ginI 
i2005l) has advocated that the faint auroral line flux can be 
computed from the fluxes in the strong nebular lines via the 
ff - relation. Then, using the obtained flux in the auroral 
line, accurate oxygen abundances can be derived using the 
classical Tc - method. We will estimate (0/H)c in the spiral 
galaxies reported to be the oxygen- richest ones, NGC 3184, 
NGC 3351, NGC 5194, and NGC 6744. For comparison, we 
will also consider NGC 2903, NGC 2997, NGC 5236 as weU 
the Milky Way Galaxy. 

We describe the method used to determine oxygen 
abundances in Hll regions in Section 2. The oxygen 
abundances in the spiral galaxies NGC 3184, NGC 3351, 
NGC 5194, NGC 6744, NGC 2903, NGC 2997, NGC 5236 
and the Milky Way Galaxy are determined in Section 3. We 
discuss the reliability of the derived abundances and sum- 
marize our conclusions in Section 4. 

For the line fluxes, we will be using the following no- 
tations throughout the paper: R2 = I[oi i] \3727+\3729/ 1 hi3, 

R-S ~ IlOIII]\i959+\S007/lHl3, R = I[0 II I]\4363/ 1 H 13 , R23 = 

R2 + R3. With these definitions, the excitation parameter 
P can then be expressed as: P = R3/(R2+R3). 
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2 ABUNDANCE DERIVATION 

2.1 Adopted equations for the Tc method 

A two-zone model for th e temperature struc ture within the 
H II region was adopted. Ilzotov et alJ ll2005l) have recently 
published a set of equations for the determination of the 
oxygen abundance in Hll regions for a five-level atom. Ac- 
cording to those authors, the electron temperature ts within 
the [O III] zone, in units of lO^K, is given by the following 
equation 

_ 1.432 

\og{Ri/R)-\ogCT ^ ' 

where 

Ct = (8.44 - 1.09 t-j, + 0.5 tl - 0.08 tl) v 



_ 1 + 0.0004 3:3 
^ 1 + 0.044 a;3 

and 

ir,-4 .-1/2 
X3. — 10 nc% 



(2) 
(3) 

(4) 



As for the ionic oxygen abundances, they are derived 
from the following equations 

12 + \0g{0+ + /H+) = log(/[oi/7]A4959+A5007//ff3) + 

1 251 

6.200 + - 0.551ogt3 -0.014t3, (5) 



12 + log(0+/jy+) = log(/[o//]A3727+A3729/-^H^) + 

gt2 - 0.034 t2 + 

log(l + 1.35a;2). (6) 



5.961 + - 0.40 log t2 - 0.034 12 + 

t2 



10 rict^ 



(7) 



Here Ue is the electron density in cm~ . 

The total oxygen abundances are then derived from the 
following equation 



o+_ 0+^ 

H+ ^ H+ 



(8) 



The electron temperature t2 of the [0 11] zone is usually 
determined from an equation which relates t2 to t^, derived 
by fitting Hll region models. Several versions of this t2 - 
ts relation have b een proposed. A widely used rel ation has 
been s uggested by ICampbell et alJ l|l986fl (see also lOarnettl 
J1992M based on the H 11 region models of lStasinskal iigS^T 
ICampbeU et alj JigSSl has found that the t2 - relationship 
can be parameterized as 



t2 = 0.7f3 + 0.3. 
It will be used here. 



(9) 



Eg. mot will be hereinafter referred to as the flux - flux or ff 
- relation. It was found that this relationship is metallicity- 
dependent at low metallicities, but becomes independent of 
metallicity (within the uncertainties of the available data) 
at metallicities higher than 12+logO/H ~ 8.25, i.e. there is 
one-to-one correspondence between the auroral and nebular 
oxygen line flux es in spectra of high-metallicity H 11 regions. 
IPilvuginl |20.Qd) derived 



log R = -4.264 + 3.087 log R2S. 



(11) 



We can now make use of the ff - relation to define a "discrep- 
ancy index", equal to the difference between the logarithm 
of the observed flux R°'"= in the [O lll]A4363 line and that of 
the flux R'''^' of that line derived from the strong [0 11] A3727, 
[O III] AA4959, 5007 lines using the ff - relation: 



Ds = log R° 



log R" 



(12) 



Since R23 — R3 /P, the ff - relation can be also expressed 
in the form R=f(R3,P). We will consider an expression of the 
type 



log _R = ai + 02 log P + 03 log R3 + 04 (log P) 



(13) 



Using a sample of H 11 regions with recent high-precision 
measurements of oxygen lines fluxes, the values of the co- 
eflicients in Eg. 1131 can be derive d. Such a sample has 
been compiled by Pilyugin (2005). iBresoUn et alJ ^00^ 
have recently detected the auroral lines [Siii]A6312 and/or 
[Nll]A5755 and determined Te-based abundances in 10 faint 
Hll regions in the spiral galaxy M 51. Unfortunately they 
were unable to detect the oxygen auror al line. Using the elec - 
tron temperature ts recommended bv lBresolin et all (|20o3) 
and the equations given above for the Te-method, we have 
estimated the value of R which co rrespond s to the ts for 
every H 11 region. It has been found jPilvuein 2005') that the 
oxygen line fluxes in six faint H 11 regions in the spiral galaxy 
M 51 (CCM 54, CCM 55, COM 57, CCM 57A, CCM 71A, 
and CCM 84A) satisfy the ff - relation. These data are in- 
cluded in the present sample to enlarge the ranges in P and 
R3. Thus, our sample consists of a total of 48 data points. 

The values of the coefficients in Eg. 11311 are derived by 
using an iteration procedure. In the first step, the relation 
is determined from all data using the least-square method. 
Then, the point with the largest deviation is rejected, and 
a new relation is derived. The iteration procedure is pur- 
sued until two successive relations have all their coefficients 
differing by less 0.001 and the absolute value of the largest 
deviation is less than 0.1 dex. The following ff - relation was 
obtained 

logP = -4.151 -3.118 log P + 2.958 log i?3 

- 0.680 (log P)^ (14) 



2.2 The ff - relation 

The fluxes R in the auroral lines are necessary to derive 
the oxygen abundances in H 11 regions using the To method. 
But they are faint and often undetecta ble in spectra of high- 
metallicity Hll regions. It was shown llPilvuginll200l that 
the flux R in the auroral line is related to the total flux R23 
in the strong nebular lines through a relation of the type 

logP = a + fexlogP23. (10) 



2.3 Characteristics of the ff — relation 

The top panel in Fig. shows the flux R in the oxygen 
auroral line as a function of the flux R3 in strong nebu- 
lar [Olll]A4959 + A5007 lines. The open and filled circles 
show individual Hll regions with 12+logO/H > 8.25. The 
relations corresponding to Eg. (1141 for different values of the 
excitation parameter are shown by the solid (P = 1.0), long- 
dashed (P = 0.7), short-dashed (P = 0.3), and dotted (P = 
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Figure 1. The ff — relations for Hll regions. Top panel. Flux R 
in the oxygen auroral line as a function of the flux R3 in the 
strong nebular [Oin]A4959 + A5007 lines for Hll regions. The 
circles (open and filled) show Hll regions with 12-|-logO/H > 
8.25. The relations corresponding to Eg. 1141 for different values 
of the excitation parameter are shown by the solid (P = 1.0), 
long-dashed (P = 0.7), short-dashed (P = 0.3), and dotted (P = 
0.1) lines. Only filled circles are used in deriving Eg. 1141 . Bot- 
tom panel. Flux R in the oxygen auroral line as a function of 
the total flux R23 in the strong nebular [On]A3727 -I- A3729 and 
[Oni]A4959 -I- A5007 lines. The open and filled circles have the 
same meaning as in the top panel. The relations corresponding to 
Eq. 1141 for different values of the excitation parameter are shown 
by the solid (P = 1.0), long-dashed (P = 0.3), short-dashed (P 
= 0.2), and dotted (P = 0.1) lines. The relation corresponding to 
Eg. 111! is shown by plus signs. 



0.1) lines. The filled circles show the data used in deriving 
Eg. 11411 (40 out of 48 original data points). 

The deviations from the S - relation given by Eg. 11411 
as parameterized by the discrepancy index Dff are shown in 
Fig. 121 as a function of oxygen abundance (top panel) and 
of excitation parameter P (bottom panel). The filled circles 
are Hll regions used in deriving Eg. 11411 . The open circles 
are the other H 11 regions. Fig. |21 shows that the deviations 
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Figure 2. Deviations from the ff — relation as parameterized by 
the discrepancy index Dg as a function of oxygen abundance (top 
panel) and of excitation parameter P (bottom panel). The filled 
circles are H n regions used in deriving the ff — relation. The open 
circles are other H II regions. 



do not show a correlation either with metallicity or with the 
excitation parameter. 

Let us compare the two ff - relations given by Eg. lilt 
and Eg. 1141 . The bottom panel in Fig.^shows the flux R in 
the oxygen auroral line as a function of the total flux R23 in 
strong nebular [O ll]A3727 + A3729 and [O ill] A4959 + A5007 
lines. The filled circles are the Hll regions used in deriving 
Eg. 1141 . The open circles are the other Hll regions. The 
relations corresponding to Eg. 1141 for different values of the 
excitation parameter are shown by the solid (P — 1.0), long- 
dashed (P = 0.3), short-dashed (P = 0.2), and dotted (P = 
0.1) lines. The relation corresponding to Eg. 1111 is shown 
by the plus signs. Inspection of the bottom panel of Fig. 
shows that the relations given by Eg. 11411 for values of the 
excitation parameter P ranging from 1 to ~ 0.3 are close 
each to other. That led IPilvueii] (|200^ to conclude that 
there appears to be no third parameter in the ff - relation. 
Indeed, examination of the bottom panel in Fig. Q shows 
that the ff - relations given by Eq.linj and Eq.jnj are very 
similar to each other and reproduce the observational data 
well at high values of the excitation parameter, P > 0.3. At 
the same time there is an appreciable divergence between 
those two relations at low values of the excitation parameter, 
P < 0.3. For definiteness's sake, we will use the ff - relation 
given by Eg. 1141 for all values of P. 

It should be noted however that the particular form 
of the analytical expression adopted for the ff - relation 
may be questioned. We have chosen a simple form. Eg. 1131 1. 
but perhaps a more complex expression may give a better 
fit to the auroral - nebular oxygen line fluxes relationship. 
Furthermore, the coefficients in the adopted expression are 
derived using calibrating Hll regions which have a discrep- 
ancy index Dff as large as 0.1 dex in absolute value. Perhaps 
smaller absolute values of Dff may result in a more precise 
relation. It can be seen that the calibration curves give a 
satisfactory fit to the observational data. At the same time, 
Fig. shows that the available measurements for calibrat- 
ing faint Hll regions (those with low R3) are very few in 
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Figure 3. The difference Alog{0/H) = Iog(0/H)T^ - log(0/H)ff 
as a function of Iog(0/H)T„ {top panel), of excitation parameter 
P {middle panel), and of discrepancy index Dff (6o44om panel). 



number. Furthermore, the R values for the faint H ll regions 
from Bresolin et al. (2004) are not measured but estimated 
from the electron temperatures ta and, consequently, they 
are not beyond question. Clearly, high-precision measure- 
ments of oxygen lines fluxes in faint H ll regions are needed 
to check the derived ff - relation. 



2.4 The (0/H)fi abundances 

Using the Te method, both (0/H)tc abundances based 
on the measured line fluxes R°'"^, and (0/H)ff abundances 
based on the line fluxes R'^^' determined from the ff - rela- 
tion, can be derived for our sample of Hll regions. Com- 
parison of their values gives us a check on the (0/H)ft 
abundances. Fig. |3] shows the difference Alog(0/H) — 
log(0/H)T, - log(0/H)fi as a function of log(0/H)Te (top 
panel), of excitation parameter P (middle panel), and of 
discrepancy index Dff (bottom panel). Inspection of Fig. |21 
shows that the differences Alog(0/H) do not show an ap- 
preciable correlation with metallicity or with excitation pa- 
rameter. There is however an anticorrelation of Alog(0/H) 
with Dff. This can be easily understood. The discrepancy 
index appears to be an indicator of the error in the auro- 
ral line R measurements. In that case, positive values of Dff 
imply overestimated R fluxes, which result in turn in over- 
estimated electron temperatures and underestimated oxy- 



n:; 

o 



Di 8.6 

o 



8.4 



1 1 1 

M51 


-i 1 1 1 — 


— 1 1 


— 1 


oO« 








o 






- 




9 








— OgC) 8 









0» 




o 




• 







1 


I.I, 


1 , 1 


1 



R (kpc) 

Figure 4. Radial distribution of oxygen abundances in the disk 
of the spiral galaxy M 51 (NGC 5194). The open circles denote 
(0/H)ff abundances. The solid line is the linear least-square best 
fit to those dat a . The filled circles show (0/H)t„ abundances from 
iBresolin et alj i2004l) . The dashed line is their linear regression 
fit. 



gen abundances, i.e. negative Alog(0/H), just the observed 
trend. In general, the (0/H)ff abundances agree satisfac- 
torily with the (0/H)tc abundances and do not show any 
systematic trend. 

In summary, the flux in the faint auroral line 
[O lll]A4363 can be estimated reasonably well from the mea- 
sured fluxes in the strong nebular lines [O ll]A3727 + A3729 
and [Olll]A4959 + A5007 via the ff - relation. Then, using 
the derived flux in the auroral line, the oxygen abundance 
can be found through the classic Te - method. 



3 OXYGEN ABUNDANCES IN SPIRAL 
GALAXIES 

Here the derived ff - relation will be applied to derive the 
oxygen abundances in a number of spiral galaxies. First, 
M 51 (NGC 5194), the oxygen-richest spiral galaxy ((0 /H)c 
= 9.54) in the sample of lVila-Costas fc Edmund^ il992l) will 
be considered. The value of (0/II)c in M 51 was derived 
recently from determinati ons of (Q/H)Tp abundances in a 
number of Hll regions bv IBresolin et aiT ^2004^ . Compari- 
son between the radial distributions of (0/H)t„ and (0/H)ff 
abundances in the disk of M 51 provides another possibility 
to test the reliability of the latter. Second, (0/H)ff abun- 
dances in NGC 3184, NGC 3351, and NGC 6744 which have 
been reported to be the most oxygen-rich spirals, will be de- 
termined. Third, we will consider three well-observed spiral 
galaxies NGC 2903, NGC 2997, and NGC 5236 for compar- 
ison. Fourth, the ff - relation will be applied to Galactic H ll 
regions to select high-precision measurements. 



3.1 Oxygen abundances in M 51 

The compilati on of published spectr a of H ll regions in M 51 
is taken from iPilvuein et al.l i2004ft. The Hll region s with 
measured R3 and R2 fluxes from ^^resolin et alJ i2004l) (their 
Tables 1 and 6) have been added. We determine (0/H)ff 
abundances for every Hll region using the derived ff - re- 
lation. The (0/H)fr abundances in the Hll regions of M 51 
are shown as a function of galactocentric distance in Fig. |1| 
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by open circles. The solid line is the linear least-square best 
fit to those data: 

12 + \og{0/H) = 8.74 (±0.03) - 0.025 (±0.004) x R. (15) 

The distance of M 51 (d = 7.64 Mpc) and t he isophotal ra- 
dius ( R25 = 5.61 arcmin) were taken from IPilvugin et alJ 
||2QQJ). T he filled circ les arc (0/H)tc abundances deter- 
mined by iBresolin et aL L20041 . The dashed line is their 
linear regression: 

12 + \og{0/H) = 8.72 (±0.09) - 0.02 (±0.01) x R. (16) 

Fig. |1| shows that the radial distribution of (0/H)ff 
abundances agrees well with that of the (0/H)tc abun- 
dances (compare also Eg. lldH and Ea. ll6ll ). This is a strong 
argument in favor of the reliability of the (0/H)ff abun- 
dances. 

Some H 11 regions show large deviations from the general 
radial trend in oxygen abundances. Some of those deviations 
can be real, but some are probably caused by uncertainties 
in the oxygen abundance determinations. It should be em- 
phasized that the uncertainties in the (0/H)ff abundances 
are not necessarily caused by unc ertainties in the line flux 
measurements. It has been noted (|Pilvuginll20'05h that the 
If - relation gives reliable results only if two conditions are 
satisfied; i) the measured fluxes reflect their relative contri- 
butions to the radiation of the whole nebula, and ii) the H 11 
region is ionization-bounded. If these two conditions are not 
met, then the derived (0/H)ff abundances may be signifi- 
cantly in error. 

3.2 Oxygen abundances in the most oxygen-rich 
spiral galaxies 

The compilation of published spectra of H 11 regions in the 
spiral gala xies NGC 3184 NGC 3351, and NGC 6744 was 
taken from iPilvuein et alj ||20^). The (0/H)fj abundances 
in the Hll regions of NGC 3184 is shown as a function of 
galactocentric distance in the top panel of Fig. |3 The solid 
line is the linear least-square best fit to those data: 

12 + \og{0/H) = 8.73 (±0.03) - 0.035 (±0.004) x R. (17) 

The same data are shown for the Hll regions of NGC 3351 
in the middle panel of Fig. |5] The solid line is the linear 
least-square best fit to those data: 

12 + \og{0/H) = 8.74 (±0.02) - 0.023 (±0.004) x R. (18) 

Similarly, the data for the Hll regions of NGC 6744 are 
shown in the bottom panel of Fig. |S] The linear least-square 
best fit to those data is given by: 

12 + \og{0/H) = 8.75 (±0.04) - 0.022 (±0.003) x R. (19) 

The distances of the galaxies and the gal actocentric dis- 
tances of the H II regions were taken from iPilvuein et alJ 

Jiool. 

Figs. Hand (see also Eas. J15ll7l I18ll9ll ) show that 
the values of (0/H)c in these spiral galaxies are ~ 8.75. 

3.3 Oxygen abundances in other spiral galaxies 

We will consider here three well-observed spiral galaxies for 
comparison. The compilation of the published spectra of H 11 
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Figure 5. Radial distributions of (0/H)ff abundances in the spi- 
ral galaxies NGC 3184, NGC 3351, and NGC 6744. The open 
circles show individual H 11 regions. The lines are the least-square 
fits to those data. 

regions in the spiral galaxie s NGC 2903. NGC 2 997, and 
NGC 5236 was ta ken from IPilvugin et al.l (|2004h . Recent 
measurements from IBresolin et a^^^2005^ have been added. 
Oxygen abundances were calculated for every Hll region 
using the derived ff - relation. 

The (0/H)ff abundances in the Hll regions of 
NGC 2903 are shown as a function of galactocentric dis- 
tance in the top panel of Fig. |S] The solid line is the linear 
least-square best fit to those data: 

12 + \og{0/H) = 8.68 (±0.02) - 0.026 (±0.003) x R. (20) 

The middle panel of Fig. 1^1 shows the same for NGC 2997. 
The linear least-square fit to those data is: 

12 + log(0/ii") = 8.59 (±0.05) - 0.021 (±0.007) x R. (21) 

The data for NGC 5236 are shown in the bottom panel of 
Fig. |S] The linear least-square fit to those data is: 

12 + log(0/i^) = 8.62 (±0.01) - 0.022 (±0.005) x R. (22) 

The d istances of the galaxies were taken from lPilvuein et alJ 

J2004f) . 

Fig. m (see also Eqs.^^^^) shows that the (0/H)c 
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Figure 6. Radial distributions of {0/H)fi abundances in the spi- 
ral galaxies NGC 2903, NGC 2997, and NGC 5236. The open 
circles show individual H II regions. The lines are the least-square 
fits to those data. 



abundances in NGC 2903, NGC 2997, and NGC 5236 are 
smaller than 8.70, i.e. lower than in NGC 3184, NGC 3351, 
NGC 5194, and NGC 6744. 

Thus, the maximum gas-phase oxygen abundance in the 
oxygen-richest spiral galaxies is 12-flog(0/H) ^ 8.75. This 
value is significantly lower (by a factor of 5 or more) than the 
previous valu e based on the one-dimensiona l empirical cali- 
brations (e.g. [y il^-Cy^^^ E(;^munds"l992l: IZaritskv et alJ 
Il994l:lvan Zee et al.lll998l: iGarnett ii200a) . 



3.4 Oxygen abundances in the Milky Way Galaxy 

We now consider the Milky Way Galaxy. A compilation 
of published spectra of Galactic Hll regions with an au- 
roral [Olll]A4363 line has been carried out by Pilyugin et 
al. (2003). That list contains 69 individual measurements 
of 13 Hll regions in the range of galactocentric distances 
from 6.6 to 14 .8 kpc. Recent measu rements of six Galactic 
H 11 regions bv lEsteban were added to that list, 

resulting in a total of 75 measurements. The Te-based oxy- 
gen abundances were estimated using the equations above. 
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Figure 7. Radial distribution of (0/H)t^ abundances in the 
disk of the Milky Way Galaxy. Top panel. The open circles show 
(0/H)t^ for all Hll regions. The solid line is the linear least- 
square fit to those data. Bottom panel. The filled circles show 
(0/H)t„ abundances for Hll regions with an absolute value of 
the discrepancy index Dg less than 0.05 dex. The solid line is the 
linear least-square fit to those data. The open square shows the 
oxygen abundance in the interstellar gas at the solar galactocen- 
tric distance. 



The derived (0/H)Te, abundances are shown as a function of 
galactocentric distance by the open circles in the top panel 
of Fig.|7| The linear least-square best fit to those data (with 
two points with a deviation more than 0.2 dex rejected) 

12 + \og{0/H) = 8.79 (±0.05) - 0.034 (±0.005) x R (23) 

is shown by the solid line in the top panel of Fig. |7| 

There is a relative large number of measurements of 
abundances in the Milky Way Galaxy. They show a large 
scatter since the quality of the data obtained over a pe- 
riod of more than thirty years is necessarily heterogeneous. 
Some measurements have a low accuracy. The ff - relation 
provides a way to s elect out the H ll regions w ith high qual- 
ity measurements JPilvuein fc Thuan I lioosl) . The bottom 
panel of Fig.|7|shows the (0/H)to abundances as a function 
of galactocentric distance for objects with an absolute value 
of the discrepancy index Dfj less than 0.05 dex. The linear 
least-square best fit to those data 

12 ± log(0/i^) = 8.70 (±0.04) - 0.022 (±0.004) x R (24) 

is shown by the solid line in the bottom panel of Fig. Q The 
open square is the oxygen abundance of the interstellar gas 
at the solar galactocentric distance (see below). 

The central oxygen abundance in the Milky Way ob- 
tained here, (0/H)c ~ 8.70, is significanly lowe r than the 
widely used (0/H)c = 9.38 of lShaver et aIllll983D . Our oxy- 
gen abun dance radial gradient is in agreement with that of 
[ Dafl on fc Cunha (2004). They derived the relation 

12 + \o^{0/H) = 8.762 (±0.105) - 0.031 (±0.012) x R (25) 
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from abundance determinations in young OB stars. 

We note tliat the central oxygen abundance obtained 
here for the Milky Way is close to that in other large spirals. 



4 DISCUSSION AND CONCLUSIONS 

The (0/H)fr abundances determined here rely on the valid- 
ity of the classic - method. This has been questioned for 
the high-metallicity regime in a number of investigations by 
comparison with Hll region photoionization models. One 
should keep in mind however that the existing numerical 
models of H ll regions are far from being perfect iStasihskal 
|2004) . and, hence, the statement that Hll region models 
provide more realistic abundances as com pared to the T e 
- metho d should not be taken for granted JPilvuginll2o'o3) . 
Recently Stasihska ( 2005) has examined the biases involved 
in Tc-based abundance determinations at high metallicities. 
She has found that, as long as the metallicity is low, the 
derived (0/H)to value is very close to the real one. Discrep- 
ancies appear around 12-|-log(0/H) = 8.6 - 8.7, and may 
become very la rge as the metallicity increases (Fig. la in 
IStasihskallioOsI) . The derived (0/H)t„ values are smaller 
than the real ones, sometimes by enormous factors. This 
means that, if such metal-rich Hll regions exist, the classic 
Te - method will always lead to systematically lower derived 
oxygen abundances. 

Can the effect discussed by Stasiiiska be responsible 
for the low (0/H)Te abundances obtained here? Let us as- 
sume that it is the case, and that the true central oxy- 
gen abundances in the spiral galaxies are higher than 12 
+ log(0/H) = 9.0. What radial distributions of (0/H)t„ 
abundances in the disk of spiral galaxies would one expect 
in this case? Starting from the periphery of a galaxy, the 
(0/H)Te abundances would increase with decreasing galac- 
tocentric distance (following the true abundances) until the 
(0/H)t„ abundance reaches the value of 12-|-log(0/H) 
8.7. After that the (0/H)to abundances would decrease with 
decreasing galactocentric d istance because of the Stasihska 
(Fig. la in lStasiiiskal l2005') effect, although the true abun- 
dances continue to increase with decreasing galactocentric 
distance. Thus, the radial distribution of (0/H)tc abun- 
dances should show a bow-shaped curve with a maximum 
value of 12-flog(0/H) ~ 8.7 at some galactocentric distance. 

Do the actual data (Fig. 5-7) show such a behavior? 
The derived radial distributions of (0/H)to abundances in 
the disks of spiral galaxies do not show an appreciable cur- 
vature, and the (0/H)tc abundances increase more or less 
monotonically with decreasing galactocentric distance. Only 
the H II regions in the very central parts of the spiral galax- 
ies NGC 3351 (the middle panel of Fig. |3 and NGC 2997 
(the middle panel of Fig.|SJ may show (0/H)to abundances 
which are slightly lower than expected from the general ra- 
dial abundance trends, but the effect is very slight and may 
not be real if errors of 0.1 dex (Fig. 3) in the calibration 
and in the line intensity measurements are taken into ac- 
count. To settle the question of whether the abundances in 
the metal- richest H ll regions are affected by the Stasiiiska's 
effect, many more accurate measurements of H ll regions (in- 
cluding those in the very central parts) in the most oxygen- 
rich spiral galaxies are necessary. 

Thus, while our results do not rule out the possible ex- 



istence of the Stasiriska's effect, they suggest that the great 
majority of H ll regions in galaxies are in a metallicity range 
where this effect is not important. Indeed, the Stasihska's 
effect becomes important only at oxygen abundances higher 
than 12 + log(0/H) ~ 8.7. We have found that this "crit- 
ical" value is reached only in the central part of the most 
oxygen-rich galaxies. We conclude that, in the light of our 
present results, Hll regions with 12-|-log(0/H) > 8.7 are 
rare, if they exist at all. It should be noted that Stasihska 
also questione d the existence o f such metal-rich H ll regions 
(Section 2.2 in IStasihskall2005ll . 

Another factor that can affect To-based abundance de- 
termin ations is the te mperature fluctuations inside Hll re- 
gions ^|Peimber1^ll967^ . If they are important, the (0/H)to 
abundance would be a lower limit. In most c ases the effe ct 
appears to be small, probably under 0.1 dex iPagelll2003l) . 

Fortunately, there is a way to verify the validity of the 
Tc - method at high metallicities (iPilvugin 2003'). High- 
resolution observations of the weak interstellar OiA1356 ab- 
sorption lines towards stars allow one to determine the in- 
terstellar oxygen abundance in the solar vicinity with a 
very high precision. It should be noted that this method 
is model-independent. These observations yield a mean in- 
terstellar oxygen abundance of 284 - 390 O atoms per 10^ H 
atoms (or 12-|-log(0/H) = 8.45 - 8.591 iMever et alJll998l: 
[Sofia & Mover 2001; C artledee et al.ll2004) . lOliveira etlo] 
(2005) have determined a mean O/H ratio of 345 ±19 O 
atoms per 10** H atoms (or 12-flog(0/H) = 8.54) for the 
Local Bubble. Thus, an oxygen abundance 12-|-log(0/H) = 
8.50 in the interstellar gas at the solar galactocentric dis- 
tance seems to be a reasonable value. The value of the oxy- 
gen abundance at the solar galactocentric distance traced by 
(0/H)Te abundances in Hll regions is then in good agree- 
ment with the oxygen abundance derived with high precision 
from the interstellar absorption lines towards stars (open 
square in the bottom panel of Fig. |7|l. This is strong evi- 
dence that the classic Te - method provides accurate oxygen 
abundances in high-metallicity Hll regions. 

The measured oxygen abundance in the solar vicinity 
provides an indirect way to check our derived central oxy- 
gen abundances. The oxygen abundance in the solar vicin- 
ity is 12-|-log(O/H)=8.50, and a gas mass fraction /j, ~ 0.20 
seems appropriate. The simple model of chemical evolution 
of galaxies predicts that a decrease of /x by 0.1 results in a 
increase of the oxygen abundance by '-^ 0.12 dex in the range 
of fi from ~ 0.75 to ~ 0.05. When all the gas in the solar 
vicinity will be converted into stars, the oxygen abundance 
will increase by 0.2 - 0.25 dex and will reach a value around 
12-|-log(0/H) = 8.70 - 8.75. This value is in excellent agree- 
ment with the central oxygen abundances obtained above 
for the Milky Way and other galaxies. There is no room for 
central oxygen abundances as large as 12+logO/H = 9.00 
or greater. 

The present study suggests that there is an upper limit 
to the oxygen abundances in spiral galaxies. If true, the ex- 
isting luminosity - metallicity relation should be revisited. 
One can expect the metallicity to increase with luminos- 
ity up to 12-f(0/H) ~ 8.75, but to remain approximately 
constant for higher metallicities, resulting in a flattening of 
the luminosity - metallicity relation. This suggestion will be 
investigated in a future paper. 

The maximum gas-phase oxygen abundance in spiral 
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galaxies is 12+log(0/H) ^ 8.75. The central oxygen abun- 
dance in the Milky Way is similar to that in other large spi- 
rals. Some fra ction of the oxygen is locked into dust grains 
in H 11 regions. ISsteban et al. Il|l998h found that the fraction 
of the dust-phase oxygen abundance in the Orion nebula is 
about 0.08 dex. Then, the true gas-l-dust maximum value 
of the oxygen abundance in H ll regions of spiral galaxies is 
12+log(0/H) ~ 8.85. This value can increase up to ~ 8.90 
- 8.95 if temperature fluctuations inside Ifll regions are im- 
portant. 
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